SNI: a new mechanism for gamma ray burst. I. Weak magnetic field 
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ABSTRACT 

We propose a new mechanism for high energy gamma ray burst in Supernova type 
I (SNI) explosion. Presupernova is assumed to be a binary system comprised of a red 
giant and a white dwarf with a wind accretion. The accretion flow terminates by the 
accretion shock in the vicinity of the white dwarf at the distance of the order of the 
accretion radius. The gas inside the accretion radius constitutes the main body of a 
target for gamma ray production. 

The supernova explosion and the shock propagation in the white dwarf result in the 
hydrodynamical acceleration of the outer layers of the star. It proceeds in two stages. 
The first stage is caused by the shock propagating in the outer layers of the star; 
the second stage is connected with the adiabatic expansion of the ejected shell into 
low-density medium around the white dwarf. The spectrum of accelerated particles is 
steep and the maximum energy does not exceed 1000 GeV. 

The gamma ray burst is produced due to interaction of the accelerated particles 
with the gas in the binary system. Most of the photons have energies about 100 
MeV, the total number of the emitted photons is between 10^^ and 10^"^. The typical 
duration of the burst is ~ 1 - 3 s for ~ 100 MeV photons and 10"^ s for ~ 1 GeV 
photons. 

The gamma ray burst might have one or two precursors. The first one is produced 
during the shock breakout when the shock approaches the star surface and crosses it. 
This burst is produced by the heated gas behind the shock; the radiation is blue-shifted 
due to the relativistic motion of the shell. The second burst might be produced under 
the appropriate choice of the parameters at the stage of the adiabatic expansion of the 
shell of the accelerated matter, when the shell becomes transparent for radiation. 
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The calculations are valid in the case of a weak magnetic field. The case of strong 
magnetic field will be considered in part II (in preparation). 

Subject headings: gamma rays: burst - stars: supernovae - acceleration of particles - 
shock waves 

1. Introduction 

It could be that there are a few different types of sources of observed gamma ray bursts 
(GRB). If GRB have cosmological origin,the repeaters are sources of a different type. GRB with 
smooth and complex time profiles can be produced by two different types of sources. Probably 
burst radiation is a typical phenomenon for many types of sources. 

Our aim is not to explain the observed GRB but to search for reliable production mechanisms, 
which can be discovered in future. In this sense our paper is similar to the one by Berezinsky 
and Prilutsky (1985) where the mechanism of GRB due to i^i? annihilation at the collapse of a 
compact object was suggested. At that time this mechanism could not explain the observed GRB, 
but it was used later for explanation of GRB in the case of their cosmological origin. 

We propose a mechanism of GRB in SNI explosion. More specifically, we study SNIa 
explosions, which as widely accepted occur in a binary system comprised of a carbon-oxigen white 
dwarf with nearly Chandrasekhax mass and a late-type secondary companion star (Wheeler and 
Harkness 1990, Khokhlov et al. 1993) . 

The explosion is caused by the thermonuclear reactions due to disc accretion ( Nomoto 1982) 
or wind accretion (Munari and Renzini 1992, Kenyon et al. 1993) of the gas onto the white dwarf. 
The shock wave propagation through the outermost layer of the white dwarf accelerates a small 
part of it to relativistic and ultrarelativistic velocities as it was first suggested by Colgate and 
Johnson 1960 . 

This acceleration is connected with the increase of the shock velocity, when the shock 
propagates in a medium with decreasing density. 

For a long time the Colgate mechanism (see Colgate (1984) for a review) was one of the 

challengers for acceleration of the observed cosmic rays. Now there are convincing arguments 
based on adiabatic energy losses and on the estimates of rate of accelerated particles, that this is 
not the case (e.g. Woosley, Taam and Weaver 1986). 

It is most probably true that the Colgate mechanism cannot provide all observed pool of cosmic 
rays, but it is also true that this mechanism does work under appropriate physical conditions. In 
particular it works in the case of shock propagation through the surface of a white dwarf. The 
interaction of accelerated particles with ambient gas results in the production of gamma radiation 
through p — p collisions, with the consequent tt'^ production and decays. 
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In this paper we investigate a gamma ray burst produced due to the interaction of accelerated 
particles with the accreting gas in the regime of wind accretion . In Section 2 the space distribution 
of the gas in accretion flow is considered. The propagation of the shock in the envelope of the white 
dwarf is studied in Section 3. In Section 4 we consider the additional acceleration of particles due 
to expansion of the shell in the low-density medium and calculate the resulting energy spectrum. 
The main gamma ray burst, produced as a result of the interaction of accelerated particles with 
the accretion flow, is described in Section 5. Discussion and conclusions are presented in Section 6. 



2. Gas distribution in presupernova binary system 

We shall concentrate mostly on the Supernova type la for which explosion occurs in the 
white dwarf. As a model for a presupernova we shall consider a binary system comprised of 
a carbon-oxygen white dwarf and a red giant. The accretion to the white dwarf results in an 
explosion which occurs when the accreting gas is accumulated on the carbon-oxygen interior of 
the compact companion and triggers the detonation with consequent mass ejection. 

We shall consider the wind regime of accretion when the white dwarf accretes the gas from 
the companion's stellar wind (Kenyon et al. 1993, Munari and Renzini 1992) . In this regime of 
accretion an explosion is possible for a carbon-oxygen white dwarf with initial mass 0.5 — 1.3 Mq 
and for accretion rate of helium gas Ma < 4 x 10~^ M^lyr (Fujimoto h Taam 1982, Woosley h 
Weaver 1993). . 

Observationally, such systems correspond to the symbiotic stars. The observed orbital periods 

for these stars (1 — 10 yr) imply a separation of the binary components d ~ (2 — 9) X 10^^ cm. 
The secondary component generally does not fill its Roche lobe and the mass loss is provided by 
the stellar wind at the of rate Myi/ ^ 10~^ M^lyr (Kenyon et al. 1993) . 

As a first approximation, the accretion from the stellar wind can be described by the classical 
theory of Hoyle and Littleton (1939). The parameters of this theory are the mass of the accreting 
star Ms, the density of the gas in the wind pw and its velocity in respect to the star uw- The 
theory operates with the accretion radius Race and the accretion rate Mhl given by: 

and: 

Mhl ='n-RlccPw{d)uw{d) (2) 

where G is the gravitational constant and d is the separation of the stars. We shall assume the 
supersonic flow of the gas at the large distance from the star r ^ d. 

The recent numerical calculations for two-dimensional axisymmetric wind accretion were 
performed by Matsuda, Inoue and Sawada (1987), Pryxell and Taam (1988), Taam and Pryxel 
(1988,1989) , Taam, Fu and Pryxell (1991). The hydrodynamical picture of the gas flow includes 
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the shock, the region of stagnation Uacc — 0, a considerable density increase in the wake and some 
other features. The flow critically depends on the accretion rate Ma- For small accretion rates 
Ma « M^dd the flow is steady. For the large it probably exibits a time variation. The details, 
essential for our calculations, are as follows. The shock is produced at distance Rsh = (.Race with 
^ < 1 and typically ^ = 0.1 — 0.2. The accretion rate is Ma = ijMhl where rj can reach the large 
fraction rj ~ 0.6. 

Inspired by these results we adopt the following simplified picture for the binary system and 
for the gas flow in it (Fig. 1). 

The compact star (white dwarf) and the secondary component (red giant) are separated by 
distance d. The red giant does not fill its Roche lobe and its mass loss is due to stellar wind with 
supersonic velocity uw- For our calculations it is enough to assume the spherically-symmetrical 
shock with radius: 

Rsh = CRacc- (3) 

Outside the shock the density distribution is: 

Mw 
AiruwrRG 

where rjic is the distance from the red giant. Inside the shock we assume: 



PwirRc) = , ^2 (^) 



Pa{r)=Pa{Rsh)ir/Rsh)-^^^ (5) 

according to the calculations of Taam, Fu and Fryxel (1991) . Eq. (|5|) corresponds to the 
quasi-spherical accretion inside the shock with infall velocity proportional to the free-fall velocity: 

Ua{r) = Cuff{r), (6) 

where Uff = (2GMs /r)^^"^ and C = r]/(16(^^^'^). In our calculations we shall use r] = 0.5, = 0.2, 
and thus C = 0.35. 

We assume the density jump on the surface of the shock as in the strong shock: 

Pa{Rsh) =^Pw{d) (7) 

The numerical values which we use in the calculations are presented in Table 1. We shall 
summarize the notation used in Table 1: Rs and Ms are the radius and the mass of the white 
dwarf, respectively, po is given by Eq. (^), d is the separation between the stars in the binary, 
Rsh is the radius of the accretion shock, is the rate of accretion onto the white dwarf, Myy 
is the rate of mass-loss by the secondary component of the binary, uw is the wind velocity and 
Xsh = p{r)dr is the column density of accreting gas inside the shock. 
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3. The relativistic shock in SNI and the precursor gamma burst 

The explosive nuclear burning in the white dwarf produces an outward-moving shock. In this 
section we shall study the propagation of the shock in the outer layers of the white dwarf and 
beyond. 

The process of thermonuclear burning propagates in the compact star either in the form 
of a supersonic detonation wave with velocity (1.1 — 1.4) x 10^ cm/s, or in the regime of a 
subsonic deflagration wave (Khokhlov, Miiller & Hoflich 1993) .The latter can be converted into 
a detonation wave in the outer layers of the star where the pressure gradient is high. We shall 
limit our consideration in the propagation of the shock in the periphery of the star, where the 
typical shock velocity is Ugh ^ 10^ cm/s. The density of the gas in the considered outer layer is 
characterized by a large gradient. In our estimates we shall use the polytropic model with n = 3 
for which the density is 

P{r) = Po(l - (8) 

where po is a constant and Rs is the radius of the star. The propagation of the shock in the 
outermost layers of the white dwarf with large density gradient results in the growth of the shock 
velocity up to the relativistic regime (Colgate and Johnson 1960) . The velocity of the strong 
adiabatic shock wave propagation in the medium with large density gradient —dln{p) / dln{r) >> 1 
in both the nonrelativistic and relativistic regimes may be approximated by the expression given 
by Gnatyk (1985) : 

^shl^sh = const X p-", (9) 

where (5sh = Ush/c is the dimensionless velocity of the shock, c is the velocity of light and 
Tsh = (1 ~ l^sh)~^^'^ Lorentz factor of the shock. The parameter a in this approximation 

is a ~ 0.2 for both nonrelativistic and relativistic shock waves. Analytical solutions for the 
nonrelativistic shock wave propagation in exponential and polytropic density distribution give the 
values 0.17 < a < 0.23 for abiabatic index 4/3 < jad ^ 5/3 and polytropic index 1.5 < n < 3.25. 
For the ultrarelativistic shock Johnson and McKee (1971) obtained analytically a = 0.232. In our 
calculation we shall consider two extreme possibilities: a = 0.2 and a = 0.232. 

Let us calculate now the maximum value of TshPsh for the shock in the star. Following 
Colgate (1984) we shall use instead of r another variable F, which is the fraction of the star mass 
beyond radius r. Using eq.(|8|) one obtains: 



P = p'^'t^s- (10) 



Then introducing the initial condition T^hPsh = (XshPsh)i at F = Fi we obtain for the shock 
velocity in the shell: 

TshPsh = {rshf3shU^)--^''. (11) 
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The numerical calculations (Khokhlov, Miiller and Hoflich 1993 and reference therein) show that 
the shock starts to be accelerated effectively at the depth Fi ~ 0.05. This is caused by the large 
density gradient in the outer layers of the white dwarf. Hereafter we shall use = 0.05 as the 
initial value. According to the above mentioned numerical calculations the velocity of the shock 
at F = Fi is 

{f3sh)i = 5 X lO-^A;,, (12) 

with ks limited between 0.7 and 1.0. We shall use for our calculations kg = 1.0. 

The maximum value of Tg^Psh corresponds to the smallest value of F at which the shock 
still exists in the outer highly rarefied layer. The shock breaks out at the radius Vmax where the 
dissipation of energy due to the escape of photons or other particles becomes essential. In the case 
of radiation-dominated nonrelativistic SNII shock, this condition can be expressed as an equality 
between the optical thickness of photons and the front thickness (Imshennik and Nadyozin 1989, 
Ensman and Burrows 1992) . In the case of relativistic SN shock, the condition for the shock 
break out can be taken as the equality of the pathlenght of the nucleus in the shock transition 
zone and the residual thickness of this layer (Colgate 1974, 1984) Using (P) and (|^) one obtains: 



where Xint is the interaction pathlength of nuclei in g/cm? . Substituting this value of Fmin into 
Eq.([Tl|) one obtains: 



3, 



(rshPsh)max = {'rshPsh)i{ TTT^ ) (1^) 

Msri 

The value of Xint depends on the physical conditions in the shock transition zone. For 
the relativistic shock waves considered here, the main source of shock dissipation is ion-lepton 
dynamical friction (Colgate 1974, Weaver 1976). In the compact white dwarf case, the temperature 
at the shock front nearby the surface is of the order of (3 — 4) x W^K (see below) and the number 
density of electron-positron pairs reaches 10^'^ — 10^^ cm~^. The main energy losses of ions in this 
case are due to elastic nucleus-electron scattering and bremstrahlung radiation. The pathlength 
of ions is 1 — 3 g/cm?. In the extended white dwarf case the temperature at the surface is lower 
and the energy losses are dominated by pion production. However, deeper in the star at a; ~ 10 
g/cm?, the temperature is again high enough for e'^e" production and the effective pathlength 
can be taken as Xint ~ 10 g/cm? . Fortunatelly according to eq. (U) {Tsh/3sh)max weakly depends 
on Xint, typically as x~^^^^. In the calculations below we use Xint = 3 g/cm^. Eqs.([Tl|) and (14) 



describe the propagation of the shock through the outer layers of the star. Namely, eq.(ll) gives 
the value of F/J for the shock at the various positions F and eq.(|l4|) gives the maximum value of 
F/3 , at the moment when the shock reaches the surface of the star. 



Now we shall proceed with the calculations of the velocity of the gas and its thermodynamical 
characteristics behind the shock. After the shock crosses a layer of the gas, this layer acquires a 



-7- 



radial velocity which is somewhat smaller than the shock velocity. The gas behind the shock is 
also heated to a high temperature. This is the first stage of hydro dynamical acceleration. 

Let us give the typical values of {TshPsh)max in our models. For the compact white dwarf 
model with Ms = IAMq, Rs = I x W^cm and Po = 1 x 10^"^ g/cm^ (Table 2), Fmin from eq.(^ 
is 4.5 X 10~^^ and hence (T shl^shjmax = 9.0 and 20 for a = 0.2 and a = 0.232, respectively. For 
the extended white dwarf model with Ms = I.OMq, Rs = 1 x 10^ cm and po = 4 x lO^g/cm^ 
(Table 1) we obtain Fmin = 6.3 x 10~^^ and {TshPsh)max = 4.3 and 8.8 for a = 0.2 and a = 0.232, 
respectively. Note that these values are considerably smaller than the ones in the calculations of 
Colgate (1984) . 

Heating and acceleration of the gas behind the strong adiabatic radiation-dominated 
(jad = 4/3) shock wave propagating in nonrelativistic cold gas can be described in terms of the 
conditions on the shock front as given by Blandford and McKee (1976) : 

£2 =r2(4r2 + 3)pic2 (15) 

P2 = (4r2 + 3)pi (16) 

r^^ = (4r2 - i)'(r2 + i)/(8r2 + lo) (i7) 

where indicies 1 and 2 refer to the parameters in front and behind the shock, respectively, 
r2 = (1 — Z?!)"^'''^ is the Lorentz factor of the fluid behind the shock, p is the density of the gas 
and e is the energy density. 

The equation of state of the gas behind the shock is given by: 

P2 = ^(£2 - P2C') (18) 



From eq. (18) it is easy to obtain the following formula for the pressure: 

P2 = ^(4r2+3)(r2-l)pic2. (19) 

In our case the main contribution to the pressure is provided by the radiation 

P2 = \aKTl (20) 
for kT2 « rrieC^, and by the radiation and the electron-positron pairs 

P2 = laKTHl + l) (21) 

for kT2 >> rrieC^, where ax is a constant of the energy density of radiation. 
Hence one obtains for kT2 » rrif-c^ and kT2 « me(? respectively: 

~ ' ^ (4r2 + 3)(r2-l)pic2)i/4 (22) 



llaK 
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and 

r2 = ((4r2 + 3)(r2-i)^)^/^ (23) 

ax 

A comment to eq. (|l7| ) is in order. This equation gives the connection between the Lorentz 
factor of the shock Tgh and the Lorentz factor r2 of the fluid behind the shock. Therefore, from 
the known shock velocity at each F, given by Eq. (0), one can reconstruct using Eq.([l7D the 
distribution of Lorentz factor of the gas over the outer layers at the moment when the shock 
reached the surface of the star. The maximum value of r2 corresponds to maximum value of Tg^ 
given by eq. ( |14[ ) . Numerically for the couipact white dwavf inodGi we obtain {^2^max 

= 6.7 and 15 

for a = 0.2 and a = 0.232, respectively. For the extended white dwarf model (r2)maa; = 3.4 and 6.6 
for a = 0.2 and a = 0.232, respectively. According to eq. (|22|), the large value of (T2)max results 
in the high temperature Tl"^ at the moment when the shock reaches the surface. 

Let us discuss now the temperature of the gas behind the shock front and the precursor 
gamma burst connected with it. When the shock reaches the surface of the star, the temperature 
T2'^^ behind it is given by Eq. p^ ) with r2 = {T2)max- The corresponding values for the compact 
model are 4.9 x 10^(7.4 x 10^)^ for a = 0.2(0.232) and for the extended model 1.6 x 10^(2.2 x 10^)^ 
for a = 0.2(0.232). 

The burst of radiation emerges from the star during the shock breakout, when radiation 
energy trapped in the shock escapes through the translucent external layer of the star. The 
shock breakout phenomenon was studied in detail for the case of SNII by Grassberg, Imshennik 
and Nadyozhin (1971), Imshennik and Nadyozhin (1988, 1989), Blinnikov and Nadyozhin (1991), 
Ensman and Burrows (1992) . The picture can be described as follows. At the time the shock 
approaches the surface it becomes thick. When the precursor of the shock reaches the surface, the 
front is spread over the distance with the optical thickness r ~ 100, according to the calculations 
of Ensman and Burrows (1992) . This moment can be taken as one when the shock breakout starts 
and the radiation emerges from the star. Following qualitatively this picture it is easy to evaluate 
the total energy of the burst in the relativistic case that we are examining. In the rest system of 
the shock, this is given by the product of 4iTRg{T2^^)^ and the duration of the burst. The latter 
can be estimated as the time, Strest ~ of the photon escape from the thickness, l^, of the 
shock front. For the compact dwarf, scattering of the photons is caused mainly by e''~e~-pairs (for 
the extended dwarf their contribution is rather small). Then for optical thickness r ~ 100 is 
given by ~ r(oTng+g-)~^ ~ 2.5 • 10^(2.5) cm for a = 0.2(0.232). Taking into account the blue 
shift of the emitted photons, the fluency can be found as 

= iT2)max^nR^MTr)%/i^^r^c) 

where r is the distance to the source. Numerically this fluency is ~ 2 x 10^® erg /cm?' 
for r = 10 kpc and the compact dwarf model with a = 0.2. For the extended white dwarf 
model this flux is ~ 10^ times larger. The typical energy of the photons in the burst is 
E"^" ~ (r2) max /cTI"'' ~ 300(1000) keV for the case of compact white dwarf with a = 0.2(0.232). 
The values of the temperature Tl"** and the mean energy of the photons from the precursor burst 
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are tabulated in Table 2. The duration of the burst, 5tobs: as observed by a distant observer is 
very short. It is determined by the time delay of the photons arriving from the different points of 
the star surface. The beaming effect makes the duration shorter in the relativistic case (see e.g. 
Rees 1968). Since the emitting angle is ~ 1/r in the laboratory frame, the time delay of the 
photons from the ^ = and 9 = 6e directions is 

(StUs ~ V2ff ^^^^ 
which numerically is 4 x 10^^(7 x 10"^) s for a = 0.2(0.232). 

The accurate calculations are needed for the detailed description. However, the gamma ray 
precursor due to the shock breakout seems to be detectable by such instruments as BATSE. 



4. Hydro dynamical acceleration 

Hydrodynamical acceleration proceeds in two stages. During the first one, described in section 
3, the shock propagating through the gas with the falling density, increases its velocity according 
to eq.(pTl) until the fraction of the shell ahead of it reaches Fmin given by eq. (^). Starting 
from this moment the shock is decaying. It leaves behind the heated shell moving outwards with 
velocity (32 and Lorentz factor r2 given by eq. (p^). The shell now expands in the low density 
gas. The thermal energy of the shell as well as the work of the pressure of the inner layers, are 
converted into kinetic energy of the shell. This is the second, adiabatic, stage of acceleration. 

When the shock reaches the surface of the star, the shell has a wide velocity distribution 
from nonrelativistic values in the inner parts of the shell to ultrarelativistic ones in the narrow 
outer part of the shell. During the adiabatic stage, the nonrelativistic velocities of gas elements 
increases as ~ 2k]\jU2-, where Uf and U2 are the final and initial velocities, respectively, and we 
take /ctv ~ 1 ( = 0.91 for n = 3 polytrope and = 4/3 (Kazhdan and Murzina 1992)). The 
initial velocity U2 is given by Eq. ([T^). In the ultrarelativistic case 

T/ = A. (25) 

The power b as obtained in the analytical solutions is 6 = 1 + \/3 = 2.73 for the plane shock 
wave (Johnson & McKee 1971) and b = 2.0 for the spherical one (Eltgroth 1972) . 

In our calculations we shall adopt 2.0 < b < 2.73. 

One can approximate both relativistic and nonrelativistic cases by the relation: 

Tf = T'2 + {4k%-b)^^ (26) 

J- 2 
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Therefore, the gas element located m the pre-shock shell at a distance r, which we describe by 
the fraction of external mass F given by eq. (p!o|), is accelerated first to the Lorentz factor r2(-F) , 
given by eqs. (|ri|) , (p!7|) and then to the Lorentz factor Tf given by eq. ([2^). 



The number of nucleons accelerated to a kinetic energy per nucleon higher than 
Ek = {T f — l)mHC^ is given by 

Q^O Ek) = (-^(-^^)--^— )^- . (27) 
ruH 

Therefore, Qn{> Ek) gives us the integral spectrum of hydrodynamically accelerated nuclei. 
For the nonrelativistic case T f — 1 « 1 we have: 

F{Ek) = F,{EK/Ei)-^'-"r (28) 



where 7^"* = 2/(3a) and numerically 7,™* = 3.3(2.9) for a = 0.2(0.232). The energy parameter Ei 
is given by: 

72 

Ei = —{Tf3)fmHc'^ = AM eV /nucleon. (29) 
49 

For ultrarelativistic case >> 1 one obtains 



F{E) = Ki^^y^rf (30) 

where K = 1.1 x 10^^^ and K = 2.3 x 10^^'^ for a = 0.2 and a = 0.232, respectively. The exponent 
of the spectrum in ultrarelativistic case is 7*"* = 4/(3a6) and numerically it is equal to 7*"* = 3.3 
for a = 0.2 and b = 2.0; and 7*7* = 2.1 for a = 0.232 and b = 2.73 

Note that in all cases the values of 7 are given for integral spectra. 

The parameters of the spectra of hydrodynamically accelerated particles are given in Table 

3. 

Table 3 shows that in the compact white dwarf case the Lorentz factor of hydrodynamically 
accelerated nuclei reaches the values of ~ 100 — 1000, while for the extended white dwarf the 
maximum values of Lorentz factors are less, Fj ~ (10 — 100). The fraction 10~^ — lO"^'^ of the 
total mass of the star is accelerated to the Lorentz factors Tj > 2. The net energy of these 
particles reaches 10^^ — 10^^ ergs. The net energy of the subrelativistic particles with Ek > 500 
MeV/ nucleon is 10^^ — 10^^ ergs. The total number of nucleons Q^r accelerated to energies higher 
than Ek is given in Table 3. The particles with these energies can generate gamma radiation 
through production of neutral pions. E^ in Table 3 refers to the end of the power low spectrum. 

During the expansion of the relativistic shell its optical depth r decreases. When r >> 1 
the energy losses of the shell due to radiation are small. The thermal energy of the shell and the 
work of the pressure of the internal layers are converted into kinetic energy of the shell and these 
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processes are taken into account in the above mentioned calculations. Due to the work of the inner 
layers, the total energy of the outer layer in the end of the adiabatic expansion is larger than its 
total energy at the moment when the shock crossed it. 

When the optical depth of the shell decreases to r < 1, the remaining thermal energy will 
be mostly radiated away and at some conditions a second precursor gamma ray burst might be 
produced. The physics of this second burst is similar to that of the prompt precursor of gamma 
burst in the fireball models (Goodman 1986; Paczynsky 1986; Shemi and Piran 1990; Meszaros 
and Rees 1993) . 



5. Gamma ray burst 

The expansion of the relativistic shell after the adiabatic phase looks like the propagation 
of a powerful beam of relativistic particles. They have energies up to 1000 GeV/nucleon and 
total energy 10^^ — 10^^ ergs. These particles interact with the accreting gas producing the pions 
and gamma rays from neutral pion decays. The distribution of the gas around the white dwarf 
is described by eq. (^ for r > Rgh and eq. (|5|) for r < Rgh- The numerical parameters for this 
distribution are listed in Table 1. 

Our calculations for production of photons (through — > vr*^ — > 77) are performed for 
two energy regions: E™"^ < Ek < Ec and Ek > Ec, where E'^^'^ = OAGeV/nudeon and 
Ec = WGeV/nucleon. 

For the low energy region 0.4 — 10 GeV/nucleon we calculate the total number of photons 
of all energies produced in the burst at angle 9 relative to the direction to the secondary 
companion: 



^ = 2 r- ,E,9if!^ < > (31) 

where n]\[{r,6) is the number density for the nucleons in the gas, Qn{Ek) is the total number 
of accelerated nucleons with energy Ex per nucleon and < £,a{EK) > is the cross section for vr'' 
production in NN— collisions. For this cross-section we use the parametrization of experimental 
data for p + p ^ tt" + all given by Dermer (1986) . 

For energies E > 10 GeV/nucleon we shall use the scaling approximation for the cross section: 



dE-,^ E. 

where x = Ej^/Ek, -Ett is the pion energy, do = 2>2mh and: 



da{Ej,,E^) (32) 
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U{x) = 0.67(1 - xf-^ + 0.5e-^^^ (33) 

This approximation describes very well the experimental data. The differential spectrum of 
produced gamma radiation is calculated as 

— = cro drnN{r,0) — / -^Uix) (34) 

dQ,dE^ Jrs Je^ Ek 47r Je^/Ek 

where is the maximum energy per nucleon in the spectrum of accelerated protons. The 

results of the calculations are displayed in Figs. 2-4. 

In Fig. 2 the column density of gas in the binary system is plotted as a function of angle 
6 relative to the line between the centers of the stars. One can see that for all angles the main 
contribution comes from the region r < Rgh (the horizontal lines). Only at small angles the 
contribution of the stellar wind is appreciable, however for ^ > 5" it does not exceed 10%. 

The column density in both our models is considerably less than the nuclear interaction length 
(xat ~ 40 — bQg/cm?) and high-energy gamma absorption length {xahs ~ 70g/crn^). Therefore, the 
column density dependence x{9) shown in the Fig. 2, coincides exactly with the total photon flux 
dependence N^{6). 

More explicitly the distribution of the gas in the binary system is displayed in Fig. 3, where 
the density of the gas, p, is plotted against the distance from the white dwarf. The position of the 
shock with the transition to the stellar wind density is clearly seen there. 

In Fig. 3 the density profile coincides with the time profile of the gamma burst. The energy- 
integrated number of produced photons dN^/dil for given 6 (angle relative to the binary axis) 
is shown as a function of the production time t in the binary system. This time is equal to the 
time-of-flight, t = {r — Rs)/vn, for the nuclei with energy per nucleon Ek at the maximum of the 
production rate < ^a{EK) > N{Ek)- 

The flux decreases with time approximatively as {t + to)~^^^, where tp = Rs/vn- Duration 
of the burst for a remote observer is different at different energies of photons. The peak in 
the spectrum, centered at E^ ~ 70 AleV is produced by pions with small kinetic energies and 
therefore photons are emitted isotropically. The observed duration is determined by difference 
in time-of-flight for photons arriving radially (9 = 0) and from the edge of the limb {9 ~ R/r), 
where R = \/RsRsh is the mean length of the pp-interaction region and r is the distance between 
the observer and the source. It gives the duration (Af)o6s ~ 1 s{Z s) for the compact (extended) 
model. For the high energy photons produced by the protons with F S> 1 the burst is considerably 
shorter. The simple kinematic estimates show that the photons are on average emitted at an angle 
~ 0.3/r relative to the direction of the proton propagation (i.e. the radial direction). Then the 
duration, determined again as the maximum difference of time-of-flight, is 
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{At)obs = -- «o.o5-r-2. 

c 2 c 

In particular for photons with ~ 1 GeV the observed duration of the burst is 
{At)obs = 5 • 10~^(1 • 10"'^) s for the compact (extended) model. 

The spectrum of the high energy gamma radiation is very steep: it has the same power- 
exponent 7*"* + 1 as the differential spectrum of accelerated particles. At ~ 70 MeV it has the 
well-known peak from the decays of pions at rest. Therefore, the gamma burst is detected most 
effectively at Ey ~ lOOMeV. However, at E^ > IGeV the signal is also detectable. The total 
number of emitted photons varies from ~ 4.5 x 10^"^ (for the compact white dwarf model) to 
Ny ~ 3.5 X 10^^ (for the extended white dwarf model) (Table 4). The corresponding fluences are in 
the range from 0.3 to 40 photons / cm? for a distance of 10 kpc. 



6. Discussion and Conclusions 

We study SNIa explosion in a binary system comprised of a compact star (white dwarf) and a 
late- type secondary companion (red giant). The accretion onto a white dwarf causes the ignition 
of nuclear reactions and the ejection of the shell. The accretion is assumed to be in the wind 
regime with the Roche lobe of the secondary companion not filled. The accretion flow terminates 
by the shock at distance r = Rsh from the white dwarf. The radius of the shock Rgh is of the 
order of the accretion radius Race, see eq. (|l|). At r < Rgh we assume a radial accretion flow to 
the white dwarf. 

Gas within the shock radius comprises the main fraction of the target for gamma ray 
production. Its density and column density are determined mostly by the mass of the compact 
star Ms and by the parameters of the stellar wind, its density pw and velocity uw- We take the 
density distribution and the radius of the shock Rgh according to the numerical calculations by 
Taam, Fu and Fryxell (1991) . The values of the column density of the gas as a function of 6, the 
angle relative to the line between the centers of the two stars, are given in Fig. 2. The mean value 
of the column density is x = 3.4 g/cm^ for the compact white dwarf model and x = 0.67 g/cm? for 
the extended white dwarf model (Table 1). 

The particles in our model are accelerated hydro dynamically. We distinguish two stages of 
acceleration. The first one is connected with the propagation of the shock through the shell. 
The acceleration occurs if the density in the shell decreases sharply outwards: e.g. p(r) ~ r^™ 
with m > 3. In this case the shock propagates with increasing velocity leaving the hot gas with 
large radial velocity behind the shock front. Following Colgate (1984) , we describe the position 
of the considered layer by the amount of the gas in front of it. Namely we use the fraction 
F = M(> r)/Ms, where M{> r) is the mass of the shell at the distance greater than r and Mg 
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is the mass of the star. The Lorentz factor of the gas in the shell at given F can be found from 
eqs. (11), ([T7|). The uncertainties for these calculations are mainly connected with the parameter 



a, which is defined by eq. (^); we expect that it varies from 0.2 to 0.232 (Gnatyk 1985) . 

The second stage of acceleration begins when the shock reaches the outer edge of the star. 
The shell adiabatically expands into a low-density medium and the pressure accelerates the shell. 
The thermal energy of the gas and the work of the pressure are converted into kinetic energy of 
the shell. This stage terminates when the shell becomes transparent and the thermal energy is 
radiated away. The Lorentz factor of the shell acquired at this stage is F j = , where r2 is the 
Lorentz factor at the end of the first stage and 6 is a parameter which varies from 6 = 2.0 for the 
spherical wave (Eltgroth 1972) to 6 = 2.73 for the plane wave (Johnson & McKee 1971). The 
value of h gives the largest uncertainties for the calculations of this stage of acceleration. 

The numerical predictions are given for two models: compact white dwarf {Ms = IMq and 
Rs = I X lO^cm) and extended white dwarf {Ms = IMq and Rs = 1 x 10^cm)(Table 1). The 
former case is more favorable for particle acceleration. For given values of the parameters a 
and b the spectrum of the accelerated nuclei has a power-law form with the integral exponent 
7*"* = 4/(3a6). In the compact white dwarf case with a = 0.232, b = 2.73 we obtain 7*"* = 2.1, the 
maximum energy of particles ii"""^ ~ liOOGeV/ nucleon, and the total energy of particles with 
Ek > 500MeV/nucleon ~ 3 x lO^^ergs. 



For the other extreme case of extended white dwarf and the values a = 0.2, b = 2.0, the 

Irl 



acceleration has a minimal efficiency: 7*"* = 3.3, E^""^ ~ 11 GeV/nucleon and the total energy is 
W ~ 2 X 10^5 ergs. 

Two precursor gamma ray bursts are expected. The first one is connected with the shock 
breakout when the shock traverses the outer thin layer transparent for radiation. The thermal 
radiation accumulated in the shock front is radiated away when the shock reaches the surface of 
the star. The thermal photons are blue-shifted due to the relativistic expansion of the shell and 
their typical energies reach E ~ F2fcT2 ~ 300(1000) keV for the a = 0.2(0.232) compact model and 
E ~ 47(120) keV for the extended one. Relativistic effects probably convert the Planck spectrum 
in the power-law spectrum (Shemi 1994). The total energy of the burst in the reference frame 
of the observer is Wburst ~ 2 x 10^^(2 x 10^'') ergs for the a = 0.2(0.232) compact model and 
Wburst ~ 6 X 10^''(3 X 10^^) ergs for the extended one. The duration of the burst for a distant 
observer is given by eq.(p^. This time is typically very short At ~ 10~^ s for a compact white 
dwarf. 

The second precursor burst might appear at the last stage of the adiabatic expansion of 
the shell of accelerated particles behind the shock, when this shell becomes transparent and its 
thermal energy is radiated away. 

The main gamma burst is produced due to the interaction of the accelerated nuclei with 
the accreting gas. The major part of the calculated flux corresponds to the interaction with the 
gas within the shock radius Rsh- Protons are produced due to the neutral pion decays. At high 
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energies the gamma ray spectrum is very steep: it is characterized by the same spectrum exponent 
of the accelerated particles (see eq. (^) and (^0|)). At low energy the spectrum has the usual 
peak at E^y ~ 70 MeV. Most of the photons have energies > 70 — 100 MeV. The maximum 
energies of photons are of order E"^ ~ 4 — 140 GeV depending on the values of the parameters a 
and b for the compact white dwarf model and ii^^ ~ 1 — 15 GeV for the extended white dwarf one. 
The total number of photons in the burst is A'^"* ~ 2 x 10^^ — 5 x 10^'^ for the compact white dwarf 
and Ni^°* ~ 4 x 10^^ — 1 x 10^"^ for the extended white dwarf. 

We want to emphasize here that in our calculations we used very conservative choise of the 
parameters deliberately trying not to overestimate the predictions. The accelerated particles 
traverse only about "i.bg/cvT? and O.GJg/cni^ of the matter in compact and extended models 
respectively. For acceleration we used very conservative assumptions too, which resulted in a steep 
spectrum and low maximum energies of the accelerated particles in comparison with the Colgate 
(1984) calculations. 

The predicted burst radiation can be easily detected if SNIa occurs in our Galaxy. For a 
distance 10 kpc from the Earth, the typical flux of gamma radiation at E-f > 100 MeV is 1 
phot/crn^, for E., > 10 GeV it is 10"^ - 10"^ phot/cm'^. 

The duration of the burst is different in different energy intervals. At E^ > 100 MeV it is 
of the order of 1 s. For the high energy photons at E^ > 10 GeV it reduces to 10^^ s. Thus we 
predict a novel type of gamma ray bursts with properties quite different from those presently 
observed. Since these bursts are rare phenomenon, their detection needs the patrol service. 

The hydrodynamical mechanism for gamma ray bursts which we put forward here for SNI can 
work for other objects too. The crucial condition for the hydrodynamical mechanism is a large 
density gradient which provides acceleration of the shock front. The accelerated particles produce 
a detectable gamma ray burst even in the case of relatively small column density of the target. 
As an example we can mention the accretion-induced collapse of a white dwarf to a neutron star 
(Woosley and Baron 1992 and references therein) . According to these calculations 7 • 10~® part of 
the whole dwarf mass acquires velocity higher than Uf = 38000 km/ s. Then maximum velocity 
of the shock is given by (T shPsh)max = 2.8(5.2) for a = 0.2(0.232) which is only a factor 1.5 
smaller than in the case of our extended model. Therefore, we expect here the gamma burst too. 
Another example can be given probably by the explosion of a neutron star below the minimum 
mass (Blinnikov et al 1990, Colpi, Shapiro and Teulkovsky 1993) . These examples illustrate the 
attractiveness of the hydrodynamical mechanism for astrophysical applications. 

B. Hnatyk is grateful to Laboratori Nazionali del Gran Sasso and in particular to the director 
Prof. P. Monacelli for hospitality. 
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Fig. 1. — A schematic model of a presupernova la binary system. The secondary companion (red 
giant) loses matter via stellar wind and the white dwarf accretes it in the regime of wind accretion. 

Fig. 2. — The total column density of gas in the binary system as a function of angle 9 relative to 
the binary axis (the line between the centers of two stars). 

Fig. 3. — The time profile of the gamma ray burst. On the left ordinate axis the energy-integrated 
number of photons emitted per Is and Isr is displayed . The emission time t shown in the lower 
absciss axis is defined as the time-of-flight t = {r — Rs)/vn for the nuclei with energy per nucleon 
Ek at the maximum of the production rate < ^a{EK) > N{Ek)- The time profile curve (full 
upper curve) coincides with the density distribution p[r) given by the values on the right ordinate 
axis and the upper absciss axis. The minimum distance corresponds to the surface of the star. The 
curves are given for the cases of compact and extended white dwarfs and for 9 = 90°. Note that the 
gamma radiation is almost spherically symmetrical (see Fig. 2). 

Fig. 4. — Energy spectra at E^ > IGeV for compact white dwarf (solid curves) and extended white 
dwarf (dashed curves). The corresponding values of the parameters a and b are shown in the Figure. 
The plotted quantities are the number of photons in the burst per IGeV. 
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Table 1: The values of the parameters for the white dwarf and the wind accretion flow. 



parameters 


Compact white dwarf 


Extended white dwarf 


Rs, cm 


1.0 • 10^ 


1.0 • 10^ 


Ms/Mq 


1.4 


1.0 


po,g/cm^ 


1.0 • 10^2 


4.0 • 10^ 


d, cm 


3.0 • 10^^ 


3.0 ■ 10^3 


Rsh, cm 


8.3 • 10^2 


5.9 • 10^2 


Ma,MQ/yr 


2.4 • 10-^ 


1.2 ■ 10-^ 


Mw,MQ/yr 


1.0 • 10"^ 


1.0 • 10-^ 


uw, cm/s 


3.0 • 10^ 


3.0 • 10^ 


x{r < Rsh),g/cm'^ 


3.4 


0.67 



Tabic 2: The values of the shock and gas parameters at the moment of shock wave break out. 

Compact white dwarf Extended white dwarf 

Parameters a = 0.2 a = 0.232 a = 0.2 a = 0.232 

Fmin iXTcP^ 4.5 • 10"^^ 6.3 • 10-^5 6.3 • 10"^^ 

{f3shTsh)max 9-0 21 4.3 8.8 

{T2)max 6.7 16 3.4 6.6 

Tp'',K 4.8-10^ 7.4-10^ 1.6-10^ 2.2-10^ 

ESf'',keV 280 1020 48 125 



-20- 



Table 3: Values of the parameters of hydrodynamically accelerated matter at the end of the 
acceleration phase. Numbers in the brackets for Qn give the powers of ten. 







Com,pact white dwarf 




Extended white dwarf 




a - 


= 0.2 


a = 


0.232 


a - 


= 0.2 


a = 


0.232 


Parameters 


6 = 2.0 


6 = 2.73 


6= 2.0 


6 = 2.73 


6 = 2.0 


6 = 2.73 


6 = 2.0 


6 = 2.73 




45 


170 


230 


1500 


12 


27 


44 


160 




42 


158 


214 


1400 


11 


25 


41 


150 




3.3 


2.4 


2.9 


2.1 


3.3 


2.4 


2.9 


2.1 


Qn{Ek> SOOMeV) 


3.9(48) 


3.9(48) 


4.0(49) 


4.0(49) 


2.8(48) 


2.8(48) 


2.9(49) 


2.9(49) 


QnIek> MQMeV) 


5.2(47) 


5.2(47) 


7.0(48) 


7.0(48) 


3.7(47) 


3.7(47) 


5.0(48) 


5.0(48) 


Qn{Ek > 9AGeV) 


1.4(43) 


1.8(44) 


8.6(44) 


7.5(45) 


1.0(43) 


1.3(44) 


6.1(44) 


5.4(45) 



Table 4: The parameters of the main gamma ray burst. The numbers in the brackets for N~f give 
the powers of ten. 

Compact white dwarf Extended white dwarf 

= 0.2 = 0.232 = 0.2 o = 0.232 

Parameters 6 = 2.0 6 = 2.73 6 = 2.0 6 = 2.73 6 = 2.0 6 = 2.73 6 = 2.0 6 = 2.73 

W^,ergs 2.6(42) 6.2(42) 3.2(43) 7.2(43) 5.6(41) 1.4(42) 7.0(42) 1.6(43) 

N^^^phot 1.5(46) 3.9(46) 2.0(47) 4.5(47) 3.5(45) 8.8(45) 4.4(46) 1.0(47) 

E:^,GeV 4.2 16 21 140 1.0 2.5 4.1 15 

N^{10GeV),ph/GeV 5.0(38) 2.0(40) 6.0(40) 1.0(42) 4.0(33) 4.0(37) 8.0(37) 2.0(41) 



